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1 Introduction 

Hundreds of papers dealing with "Lindblad equations" have been written since 
their fundamental importance in the theory of open quantum systems was shown 
(see |L76t IGKS76] ). They give the proper mathematical form to Markovian 
semigroups of completely positive norm continuous maps which are needed in 
this context |D76j . The early investigations (see, for example |S76t IF78t IS80j ) 
mostly aimed at establishing theories for the approach to thermal equilibrium, 
and all the subsequent investigations pursued special physical questions; none of 
them, however, dealt with the general mathematical structure. And this is now 
the theme of this paper. We give an analysis of the mathematical properties. It 
is analogous to the functional analysis of unitary groups of evolution, which was 
presented in the early time of quantum mechanics. We continue the study on 
semigroups of completely positive maps, acting on states for systems with finite 
dimensional Hilbert spaces, which we began in |BNT08] . 

While Schrodinger equations are studied in Hilbert space spanned by the 
pure-state vectors, the Lindblad equations will be studied on two levels: One is 
the linear space of density matrices, the other is the underlying Hilbert space. 
A relation between these two levels is established. In this relation we find a 
clear distinction between different processes: Decay, Dissipation and Dephasing. 
Characterizing these processes involves a structuring of the underlying Hilbert 
space H. Each Lindblad generator is related to a special decomposition of TC 
into mutually orthogonal subspaces. This decomposition can be seen as a gener- 
alization of the spectral decomposition of Schrodinger operators. Orthogonality 
of the subspaces holds in spite of non-hermiticity of the defining operators. 

We do not discuss "decoherence" . The meaning of this concept in physics is 
connected with the transition from quantumness to classicality, |J03t IZ03at IZ03b] , 
the mathematics of it requires the definition of a preferred basis. "Dephasing" 
is a similar concept, and we define what it means precisely in a mathematical 
fashion. 

Our starting point, the connection with the earlier studies, is the result of 
|GKS76l[I7f6] : 

1 PROPOSITION. Generators of semigroups: Every generator of a semi- 
group of completely positive trace preserving maps T* : p{s) t— > p{s + 1) for t > 
on the set of finite dimensional density matrices p, can be written in the form 

p = V{p) = -i[H,p\ + Y,'^uAp) (1) 

a 

where H = W is a hamiltonian. Transition operators ha define the irre- 
versible parts, the simple generators: 



Vh{p) = hph^ - \{h^hp + ph^h) 



(2) 



SGL2 BN June 18, 2008 



3 



This representation of V is known as "diagonal". For our purposes it is 
optimal. All what is needed to characterize the processes can be extracted from 
the set of the operators {H, ha, h]^}, and from the algebra of operators which 
commute with each element of this set. (A remark on the notation: We use 
here ha with lower index, instead of the use of an upper index, which has been 
employed in our first paper [BNTOS] .) It is well known, |AFni[ IBPn2j . that the 
division of "D into a sum of several simple generators, attribution of {H, ha} to 
"D, is not unique. Different sets of operators ha can be attributed to a given 
generator T>, some actions can be shifted from the Hamiltonian H to the transfer 
operators ha or in the other way. There exist attempts to favor special ways, like 
demanding Tr[/iQ,] = 0, but we make here, in this paper, no restriction. 

The mathematical analysis is of course related to special questions concerning 
physics: 

• Existence and characterization of stationary states. 

• Geometry of the paths and characterizing Hilbert space subspaces. 

• Aspects of symmetries. 

• Perturbation. 

Identifying stationary states is closely related to identifying special subspaces 
of the Hilbert space. The summary of our investigations on the first two of these 
themes can be stated as follows: 

2 THEOREM. Structuring of the Hilbert space: 

1. Decay: The Hilbert space can be represented in a unique way as a direct 
sura of two orthogonal subspaces, H = PqH © PoT^, where P^Ti is the 
maximal decaying subspace, i.e. 

Vp G S : hm P^T\p)P^ = 0, (3) 

and PqH contains no decaying sub-subspace. 

yQ<Po3peS: limsup Tr [gr* (p)g] ^ 0. (4) 

2. Dephasing: The collecting subspace PqH can further be divided in a 
unique way by splitting Pq = '^f.Qo,k, Qo,k ■ Qo,e = SkeQo,k, into minimal 
subspaces Qo^k'^ with relative dephasing 

Vp(0), VA;,£, k^l: lim Qo,fcP(t)Qo,^ = 0, (5) 

t— >oo 

and the property, that the time evolution of each block Qo,kp{t)Qo,i is inde- 
pendent of the other blocks. 
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3. Asymptotics: Each subspace Qo,k'H can be represented as Qo,k'H — 

CMk) (^T-CoQ^k; so that the time evolution at large times is described in block 
form, with Hamiltonians Hq ^ ® 1 acting on Qq^uH and with unique density 
matrices acting on pk on 7Yoo,A; inside each block: 

Vp(0) 3 {Xk, Rk} : lim \p{t) - Afee-^^°''=*i?fee^^°.'=* <8) = 0, (6) 

t—^0O 

k 

where each is a positive matrix with trace one acting on C"'*^'^^, and 

o<\k<i, EfcAfc = i. 

The set of stationary states is given by the set of density matrices 

0Afei?fe0pfe, (7) 

k 

where [i^o.fc; Rk] = 0, and pk again unique for each k. 

4- Dissipation: The limiting density matrices inside the minimal blocks are 
of maximal rank: rank(pfe) — dim(?ioo,it)- 

5. Cascades with basins: The decaying subspace PqIH. can further be 
divided by splitting it as a cascade with "basins" Pk,iH, k > 1, all mutually 
orthogonal, Pk,e-Pj,m = ^kj^mfPk/, and arranged in levels PkH = @pPk/'H. 
So Pq = Efc>i£-Pfc,£; and the time evolution is like a "flow", out of each 
basin into the "lower" levels including collecting basins Po^k'H G PqTC, where 
each Qo/Ti- mentioned in the item above contains n{€) basins: 

for t>Q: Pyr\Pk,epPk,e)-Pj 7^ if J < K (8) 

Pj-r'{Pk,ipPk,e)-Pj = if 3>k, (9) 
Pk,m-r\Pk,epPk,e)-Pk,m = if m^t (10) 

In special cases there is a possibility of unitarily reshuffling some basins, 
defining other basins Pk,i — U ■ Pk,e ■ ■ The number of minimal basins is 
unique, as are their dimensions. 

All projectors Pk, Pk,e o^nd Qo^k are orthogonal projectors. 

Perpendicular to the structuring into levels there are dissections into "enclo- 
sures" with a division of density matrices into blocks with mutually independent 
evolutions. 

3 THEOREM. Enclosures and blocks: // there exists a set of mutually 

orthogonal projectors Qm, each commuting with H and every ha, then the basins 
Pk/Ti. can be chosen in such a way that each subspace Qm'H is a direct sum of 
basins. We call such a Qm'H an enclosure. The time evolution of any density 
matrix p splits into mutually independent evolutions of blocks QmPQe- 



r\QmPQl) = QmT\p)Qt 



(11) 
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As an example for the structuring, already well known, may serve the Gro- 
trian diagram for Helium. The dissection into enclosures is there the distinction 
between ortho- and para-helium. The lowest energy levels represent the collect- 
ing subspace. Energy levels with angular momentum zero are basins, those with 
higher angular momenta J can further be split into basins, by diagonalizing a 
component of J. The choice of the axis of the component is not unique, but 
different kinds of splitting, with different basins, are unitarily equivalent. 



The Qm are projectors onto the "enclosures" - for details see Section 3.1 In 



the Heisenberg picture, with time evolution of the observables, they are invariant. 
This is a conventional symmetry. But there are cases when other symmetries turn 
up in non conventional way, without conserved observables; cases when conser- 
vation of observables appears, without an overall symmetry; and the conserved 
observables need not commute with all the ha, and they need not form an alge- 
bra. "Degeneracy", the existence of several stationary states, may be connected 
with occurrence of symmetries or not. The well known connection of symmetries 
with conservation laws, fundamental for Langrangian mechanics and quantum 
mechanics, is here no longer valid. All this is discussed in Section |4} It is also an 
important aspect when perturbations are studied, as is done in Section [6] 



2 Basic properties of the superoperators 

2.1 The Hilbert Schmidt space of matrices 

Consider the set S of states as the set of n x n density matrices pij = {i\p\j), 
with n = dim(7i), {\i)} some basis of TC, and consider it as embedded into the 
linear space of complex n x n matrices. Now consider the general mathematical 
wisdom on linear differential equations with constant coefficients, acting in a finite 
dimensional space, which is here the space spanned hj nxn dimensional matrices: 
To each superoperator V there exist eigenmatrices a, proper or generalized: 

V{a) = Xa or - A)"((t) = 0, 

So one has the special time evolutions: a{t) = e^* ■ Polynomial{t,T>"'{a)). The 
general solution to d = T'(cr) is a linear combination of these special solutions. 
But the eigenmatrices to our T> may be not self-adjoint, although V preserves 
self-adjointness. Instead, there may be pairs of eigenvalues {A, A*} with pairs 
of eigenmatrices {cr, a^}. That's the reason why we do not stay in the space of 
self-adjoint matrices. 

1 

We endow this linear space with the Hilbert Schmidt norm ||cr|| = (Tr[o"V])2 
and the inner product ((o"|p)) = Tr[(T'''p]- The HS-superspace of operators enables 
the definition of the adjoint of a super-operator. The adjoint super-operator V"^ 
generates the evolution of observables in the Heisenberg-picture, F = T>'^{F). It 
acts as 
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V\F)=z[H,F] + J2l^liF), (12) 

a 

where 

Vl{F) = h^Fh-\{h^hF + Fh)h). (13) 

We remark that self-adjointness of T> is rather exceptional. See Section [5] for 
examples. 

2.2 Evolutions inside the set of states 

T* preserves self-adjointness, trace-norm and, for t > 0, positivity; it maps S i— >■ 
S, the set of normed positive matrices into itself. 

The dissipative character of the evolution is reflected in the general changes 
of eigenvalues r(t) of p{t). Zero eigenvalues may become positive, positive eigen- 
values may increase or decrease. But the decrease is bounded from below: 

4 PROPOSITION. A differential inequality: Ifr{t) is a non-negative eigen- 
value of p{t), its change in time is bounded from below as 



r{t). (14) 



Proof. Differentiating the eigenvalue equation p(t)tp{t) = r{t)tp(t), and using 
{'ipl'ip) = gives 

a 

a 

= -5^(V^|/il/i.|^)r(t). (15) 

a 

Here, the finite dimension of the Hilbert space is essential: it implies that the 
number of simple generators which are necessary to form some given V is bounded 
by dim^(?-^) — 1. All the sums over a are finite. □ 

An immediate consequence is that the positive eigenvalues rj of the density 
operator obey the inequality 



rj{t) > exp ^- \\haft^ ■ rj(0) > 0. 



The rank of p cannot decrease in finite time, purification can only occur in the 
limit t — s> oo. 

Moreover, we can affirm quite generally 
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5 PROPOSITION. Existence of stationary states. For each V there exists 
at least one stationary state p G S. 

Proof. The density matrices in B{T-C) form a compact convex set S. We know tliat 
the semigroup maps S into itself, and we may consider the map p(0) ^ p(At) for 
some fixed time interval At. By the Brouwer fixed point theorem, there exists at 
least one fixed point, p(At) = p(0). So there exists the stationary state 



p{t)dt, (16) 



satisfying V{p) = 0. □ 



This abstract statement is needed in 3^ In the following it is then elucidated 
by more concrete formulas. 



3 Processes and Structuring 
3.1 Definitions 

Enclosure is not a standard concept, although it appears in Hamiltonian dy- 
namics as a consequence of the conservation of energy. We define "enclosure" as 
a subspace /C with the property that for all p(0) the expectation of the orthogonal 
projector onto this subspace, Tr[Picp{t)], is constant in time. In the Schrodinger 
picture this means that the system can neither leave nor enter this subspace. In 
the Heisenberg picture it means invariance of the observable Pic in the course of 
time. Each P/c is one of the projectors Qm appearing in Theorem [3j So, enclo- 
sure is a case of "conservation", but in the special way that it is an orthogonal 
projector which is the conserved observable. In Section [5] it is shown that there 
are cases where invariant observables exist, but no invariant projector, different 
to the evolutions under Hamiltonian dynamics. 

Our analysis of general Lindblad equations T* starts with investigating an 
eventual enclosing of the system in the subspaces )C with conserved projectors 
P/c- Then, inside the enclosures, decay and/or dissipation occurs. The phase- 
relations between the enclosures may eventually show a dephasing. Dephasing is 
a typical quantum-effect, but it is also related to decoherence, the transition from 
quantum to classical appearance, |Z03at [Z03bj . We introduce precise mathemat- 
ical definitions of these physical events of irreversible processes: 

6 DEFINITION. Irreversible events: 

Decay occurs, if there exists a decaying subspace fC; i.e. 



Vp(0) : Tr[P^p(t)] 0. 



(17) 
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Dissipation occurs, if there exists a subspace K, such that 

Vp(0) = Pycp(O) : rank(p(t)) ^ dim{lC). (18) 

Dephasing occurs for phase relations, in the case of conservations of Pjc and 
Pc, with K, ± C, if 

Vp(0): P,cp{t)Pc^O. (19) 
All arrows are meant as limit oft to infinity. 

Enclosures, decaying and - their complement - collecting subspaces are the 
analogue to the energy eigenspaces in Hamiltonian dynamics. 

3.2 Enclosures and conservation of projectors 

In the Heisenberg picture "conservation" of P means X>^(P) = 0. In the Schrodinger 
picture this means that the state of the system can neither leave nor enter the 
subspace PH. 

7 LEMMA. Conservation of P — P^ — P^ is possible if and only if this observ- 
able commutes both with the Hamiltonian H and with all the transfer operators 
ha', in other words, iff it is an element of the von Neumann algebra {H, ha}' . 

Proof. We consider p — PpP, systems in the subspace PTC. Leaving this subspace 
is possible for the system, if 

P^V{p)P^ = ^ P^ha p hiP^ ^0 ^ 3a: P^h^P ^ 0. 

a 

Entering the same subspace means leaving the orthogonal complement and is thus 
possible for the system if 3q; : Ph^P'^ 7^ 0. That neither leaving nor entering is 
possible implies therefore Va : Pha = PhaP = haP. Under this restriction one 
gets for the evolution of the phase relations, the off diagonal part, P pit) P^: 

PV{p)P^ - -iP[H,p]P^ = iPpPHP^, 

which vanishes for every p = PpP iff H commutes with P. And its vanishing is 
necessary, otherwise P^V^{p)P^ = P^HPpPHP^ > 0. □ 

Assume that both Pi and Pj are conserved projectors. The commutations 
imply V{Pip) = PiD{p) and V{pPj) = V{p)Pj. So the evolution of each block 
PipPj is independent from the other parts of p. The set of conserved projectors 
generates the algebra = {H, ha}'. A maximal abelian subalgebra of the von 
Neumann algebra J\f gives a set of mutually orthogonal minimal conserved pro- 
jectors Pi, a decomposition of Ti, and a decomposition of p as a block- matrix, 
with split evolutions of the blocks. Such a decomposition is unique if J\f is an 
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abelian algebra. If this algebra is not abelian, then different decompositions are 
possible. They are related by unitary transformations with U E M which leave V 
invariant. This follows from the classification of finite dimensional von Neumann 
algebras, stated for example in (JOSb j. 

Conserved projectors come along with an algebra of conserved observables, 
{iJ, ha}' 1 and, moreover, with a dynamical symmetry: 

T\U^ pU) = p{t)U if Ue{H,K}'. (20) 

But there are cases of conserved observables which do not form an algebra; and, 
on the other hand, conservation of projectors is not necessary for occurrence of 
a dynamical symmetry. This is discussed in Section |4j 

A maximal set of mutually orthogonal minimal enclosing subspaces enables 
now the discussion of processes inside. 



3.3 Collection into subspaces, dissipation inside 

In studying the time evolution in |BNT08j we identified the condition "laziness" . 
A subspace PTi is "lazy" if there is no flow out of it in first order of time. This 
condition appears here as the first part for characterizing collecting spaces, the 
orthogonal complements of decaying spaces. 

8 DEFINITION. Lazy subspaces: P7i is a lazy subspace if 

Vp = PpP : |Tr[Pp(t)P] = at t = 0. (21) 

9 LEMMA. PH is a lazy subspace 

Wa: K- P = P -K- P. (22) 

Proof. Since T* conserves the trace, one has Tr [Pp(t)P] = — Tr[P-'-p(t)P-'-]. At 
time t = this is to be expressed with the generating operator as 

-Tr[P^P(p)P^] = -Y, Tr[P^/i„p/it P^] = - ^ Tr[(P^ ■ K ■ P)p{P ■ hi ■ P^)]. 

a a 

So, if P-*- ■ ha ■ P does not vanish for each a, there exists a state p = PpP, such 
that Tr[Pp(t)P] ^ 0. This inequality holds in particular for any p = PpP which 
has rank(p) = Tr[P]. □ 

10 DEFINITION. Collecting subspaces 

PH is a collecting subspace, if 



Vt > Vp : T\PpP) = PT\PpP)P. 



(23) 
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11 LEMMA. PTi. is a collecting subspace, iff it is a lazy subspace also satisfying 

P{zH-\Y,hiK)P^ = Q. (24) 

a 

Proof. Since the Lindblad equation is of first order in taking the time-derivative, 
the equation 

Vp : V{PpP) = PV{PpP)P (25) 



is equivalent to ( 23 ) . The laziness condition is equivalent to disabling, as is shown 



above, appearance of diagonal blocks P pP : 

Vp : P^V{PpP)P^ = 0, 



and (24) is equivalent to disabling appearance of off-diagonal blocks of p{t), 

Vp : PV{PpP)P^ = 0. 
This is seen just by inserting the Lindblad equation ([T|: 

PV{PpP)P^ = PpP{iH - \ hiK)P^ (26) 

□ 

For the system there is no way (back) out of PTi. But some part of p may 
enter this subspace, if P does not commute with the Hamiltonian H and all the 
transition operators ha- In general P-T^{P^pP^)-P 7^ 0. If P does not commute 
with all the h^, this occurs for some p, for example p = P-'"/Tr[P"'"], in first order 
in time, i.e. P ■ V{P^pP^) ■ P 7^ 0. If P does commute with all the ha, but 
does not commute with the Hamiltonian, entering the subspace occurs in second 
order: P ■ {P^ pP^) ■ P 7^ 0; see the proof of Lemma [Tj 

If P^ projects onto a decaying subspace, then PTi. is collecting. This can be 



seen in the structuring of 7i, performed in Section |3.4[ For the simple V = T>h 
there was no reason for such a detailed investigation: A collecting subspace is 
there just a proper eigenspace for the eigenvalue zero of h. But here the struc- 
tures are in general richer. There is still an evolution going on in the collecting 
subspace. It can be seen as an evolution for density matrices defined on the 
reduced Hilbert space PH.. 

p = PpP =^ p = V{p) defined with = PHP, ha = Ph^P}. (27) 

If the collecting subspace is minimal, i.e. does not contain a smaller collecting 
subspace, the evolution inside is of the same type as in a minimal enclosure; it is 
dissipative, unless dim(P?i) = 1. 

The relations between characteristic subspaces and stationary states can now 
be analyzed, first for one direction of implications: 
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12 LEMMA. If p is a stationary state and P is the projector onto its range, 
then PTi is a collecting suhspace or an enclosure. 



Proof. This proof is done by revisiting the proofs for Lemma [9] and Lemma 11 



with the special extra situation p = PpP and range (p) = PH. □ 

The other relations state that inside a minimal enclosure, with no decay inside, 
there holds 

13 THEOREM. Uniqueness of minimal stationary states. Let JC be a 

suhspace which is a minimal enclosure or a minimal collecting suhspace, con- 
taining no smaller enclosure or collecting suhspace. Then there exists one and 
only one stationary state supported hy /C. Its density matrix has maximal rank, 
rank(p) = dim{)C). 

Proof There has to exist a stationary state p, see Theorem |5| applied to the 
restricted evolution inside of JC. Linearity of the Lindblad equation implies that 
the stationary states, if more than one, form a line or a (hyper)plane. This would 
include elements in the boundary of S with lower rank; but there are no invariant 
states at the houndary of S(/C): Any invariant state has a collecting subspace as 



its range (Lemma 12), and it is assumed that /C contains no smaller collecting 



subspace. □ 



3.4 Cascades of decay 

In the cases where decay occurs it might be helpful or necessary to define a struc- 
ture of the Hilbert space, analogous to the "energy levels" in Grotrian diagrams 
of atoms. 

We define the lowest level of the Cascade, PqH, as the smallest subspace 
of H. which contains all minimal collecting subspaces. The strategy of further 
procedure to give more structuring is as follows: Consider the complement of the 
lowest level, with the same evolution, except the flow out into the lowest level. 
This evolution, acting on P^^H, is generated by Vi, formed with P^HP^ and 
{PqUciPq}. Then the first higher level PiTi of the Cascade for T) is defined 
as the lowest level of the Cascade for Vi. Iteration gives a series of levels PoH, 
PiH, . . . PnH, until • PiH = H. Basins PijH in the level PiH are defined as 
minimal collecting subspaces for the reduced evolutions with no flow out into the 
lower levels. 

The precise details: 

14 PROPOSITION. Decomposition into basins. Each level PiH can he 
decomposed into a direct sum of mutually orthogonal hasins PijH. This decompo- 
sition is either unique or unique up to some unitary equivalence, which reshuffles 
basins among a set of partners with equivalent dissipations inside. 
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Proof. Construct the decomposition inductively and begin with the lowest level. 
Consider a sum of minimal collecting mutually orthogonal subspaces P7i : = 
0^ PajH, and consider another minimal collecting subspace PaTi-, not contained 



in PTi,. Pali, contains a unique stationary state a (Theorem 13). If Pa'H is 
orthogonal to PTC, define -Poj+i = Pa and proceed inductively. If PaTC is neither 
orthogonal to PTi. nor contained in it, consider the Lindblad equation acting 
on a, with restriction to the subspace spanned by adding PTi. and P^TC, in the 
block-matrix representation as stated in the Appendix. There R = PaP ^ 0, 
Q = PaP^ 7^ 0, since P^TC is not subspace of PTC but not orthogonal to it. Also 
S = P^aP^ 7^ 0, by positivity of a, and rank(5') = rank(P + Po-) — rank(P), 
since PaTC is minimal. The condition "laziness" on PTC, as stated in Lemma [91 



is Ca = P haP = 0. Together with the condition stated in equation (24), which 
is iPHP^ — ^ ^iiBa = (here with the notations as defined in Appendix, 
Section |8]), one gets for d = the part 

^ = -1 Y^^BiBo^S + SBiB^) = 0. 

a 

This implies Va : B^ = PhaP^ = 0, where is the restricted transfer 
operator, restricted to the subspace PTL + PaTC. Together with the "laziness" 
condition this means commutation of every restricted ha with P, and this implies 



moreover, again using equation (24), that also the restricted H has to commute 
with P. It follows that 

Poj+i7^ := {P^TC ® PaTC) e {PTC ® PaTC) 

is another collecting subspace, orthogonal to PTC. 

Iterating this procedure, until no other collecting subspace, no other station- 
ary state not already contained in PqTC := PTC is left, gives a decomposition of 
the lowest level. 

Now we investigate the restricted evolution of the complement of PqTC. Con- 
sider the subspace PqTL and the processes generated by {PQHPQ,PQhaPQ}. 
The collecting subspaces of this evolution, with all outflow into the lowest level 
disabled, give the basins of the first higher level PiTC = 0^ -Pi.j^- Then one 
iterates the disabling of the outflow of the remaining subspace, constructing the 
higher levels with basins PijTC until nothing more of TC remains. 

The lowest level is the subspace spanned by all possible stationary states. 
So its definition is unique. By induction, the entire decomposition into levels is 
unique. 

Now consider two different decompositions of a level, w.l.o.g. of the lowest 
level. This gives exactly the situation treated above, with at least one stationary 
state a with support in a basin PaTC which is neither orthogonal to some basin 
PijTC nor contained in it. This is the case, iff there are stationary phase relations 
PijO'Pk/- The discussion of such cases is postponed to the following Section 3.5 



using the Proposition 16 □ 
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The characterizations of basins involves a common Schur triangulation ( |L69] ) 
in block form of all the transfer operators. The blocks PijhaPk/ with i > k are 
zero. The change in time of a basin's content PijpPij consists of 

• inflow from (several) Pk,tH with k > i, 
generated by {PijHPk^i + Pk/HPij, PijhaPk/} 

• outflow into (several) Pk,iH with k < i, 
generated by {PijHPk/ + Pk/HPij, Pk/KPi^j} 

• dissipation inside the basin, 
generated by {PijHPij, PijhaPij}. 

If higher basins Pk/'H are empty, all contents of the basin Pij'H will decay if 
i>l. Dissipation leads to density matrices with full rank inside the basin; then, 
if z 7^ 0, there is some outflow since '}2k<i,i ^k/ ■ {iH + \ h\ho) ■ Pij ^ 0. 

15 PROPOSITION. Emptying of higher levels. Only the lowest level 
in the cascade carries stationary states, formed by combinations of unique states 
inside each single basin; eventually there may be stationary phase relations. The 
collection of the higher levels in the cascade, 

is the maximal decaying subspace. It is completely emptied in the course of the 
evolution. 

Proof. There are no stationary states in /C, by construction. So the diagonal 
blocks PijpPij with i > 1 vanish. By preservation of positivity, the off-diagonal 
blocks PijpPk/ and Pk/pPij have to vanish also. □ 



3.5 Dephasing and the geometry of paths 

Here we study, as t — oo, the phase relations between basins. Knowing that 
basins in the upper levels of the Cascade get empty, and phase relations involv- 
ing one or two of the decaying basins have to vanish because of preservation of 
positivity, it remains to study phase relations between minimal collecting sub- 
spaces. We may restrict the system and consider only the lowest level, PqH, 
which is collecting. In this level the time evolution is identical to the evolution 



defined in equation (27), generated by the restricted operators. Every collecting 
subspace is there, in this restriction, an enclosure. So we simplify the discus- 
sion and consider a system with can be decomposed into minimal enclosures PjTi 
without decay. The time evolution of each block PipPj is independent of all the 
other blocks. When considering the phase relations between diagonal blocks we 
may therefore simplify further, and restrict the system to a space with just two 
basins H = P/H © PjH. 
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16 PROPOSITION. Stationary phase relations. A stationary phase rela- 
tion PipPj exists if and only if there exists a unitary operator U commuting with 
H and each which intertwines between the two enclosures. 



UPj = PiU, [/2 = 1 (28) 
The stationary phase relation is unique up to a constant factor. 

Proof. Assume that such an intertwiner U exists. It creates a dynamical symme- 
try, 'D{U pU'^) = UV{p)W , and it transforms the stationary density matrices of 
the enclosures into each other, 

UpjU^ = Pi. (29) 

The commutation with the generating operators implies stationarity of phase 
relation blocks PipPj = rijUpj and PjpPi = rj^iPjU. 

On the other hand, assume that some stationary block PipPj exists. Then 
also PjpPi = {PipPjY is stationary. Since the matrices for p in diagonal blocks 
are of full rank (Theorem [Ts]), the density matrices 



(TX = UP^ + >^P^pPi + PjpPi + Pj) 

are positive for |A| small. Now at some critical value of |A| the state a\ is at the 
boundary of S, and rank((TA) < rank(Pj) + rank(Pj). This implies that Px'H, the 



range of cta, is an enclosure (Lemma 12), and that Pa £ {H, ha}' ■ Since PiTi and 



PjT^ are undecomposable and Pa does neither commute with Pj nor with P,-, the 
only possibility for such a situation is, that the von Neumann algebra {H^ha}' 
is not abelian, and 

{H, ho}' = P{H, 



including a unitary intertwiner U acting as in equation (29). □ 



Proof. Uniqueness of decomposition into basins. Proposition 16 gives 



the completion for proving Proposition 14 The decomposition into basins is 



not unique, iff there exist stationary phase relations. And such stationary phase 



relations exist, iff there is a unitary equivalence as stated in Proposition 16 □ 



The geometry of paths {pit)} C S is related to the eigenvalues of the super- 
operator v. Each path can be decomposed into a sum of at least one stationary 
state and paths of self- adjoint matrices which are eigenmatrices or pairwise sums 
of eigenmatrices of T>. 

Special paths for self-adjoint matrices: 

• eigenvalue zero <^=^ stationary state 

• imaginary eigenvalue -v^ circular path 
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• negative eigenvalue path leading straight to zero 

• complex pairs of eigenvalues with negative real part -v^^ paths formed as 
■j(j(t) + 7*0"^ (t), spiraling in to zero 

General paths arise as superpositions of special paths. Geometric considerations 
give some implications for the eigenvalues of V: Decay processes lead straight 
to the boundary of S, so negative eigenvalues are involved. Eigenmatrices whose 
range contains some part of Pfj~Ti. belong to eigenvalues with negative real part. 
The eigenmatrices whose support is in PoTi. can be chosen such that their supports 
are in single collecting basins or in blocks giving phase relations between pairs of 
basins. In this way one gets a complete set of eigenmatrices spanning the space 
of all matrices with support in PqTI. Since no path can leave S, there are no 
eigenvalues with positive real part. Moreover we can state 

17 LEMMA. To the eigenvalue zero ofV there exist only proper eigenmatrices. 
The corresponding eigenspace is spanned by positive density matrices. 

Proof. Consider the stationary eigenmatrix a. If it is not self adjoint, then also cr^ 
is an eigenmatrix, as are the selfadjoint a + a^ and ia — ia^ . If a is selfadjoint but 
not positive, consider it split as cr = p+ — p_, both parts being positive. Since T* 
is positivity preserving, both parts separately must be stationary. Now assume 
the existence of a generalized eigenmatrix Tq, with 'P(to) = cr, T>{a) = 0. Again 
one can conclude that the analog equations should hold for the adjoint matrices 
and for their linear combinations. One can therefore assume tq = t^. Integrating 
the evolution equation, assuming t(0) = tq, gives T{t) — tq + ta. Multiplying 
with some small e and adding some positive p with full rank would give a path 
starting inside S but leaving it as t gets large. This is a contradiction to the 
preservation of positivity, so no such tq can exist. □ 

The eigenvalues of V lying on the imaginary axis correspond to circular paths. 
Such cases can appear for phase relations, for off-diagonal blocks of p. 

18 THEOREM. Dephasing and non-dephasing; eigenvalues of T> 

1. For a minimal block at the diagonal, belonging to a collecting basin, {a — 
PojCT-foj}; there exists exactly one eigenmatrix to the eigenvalue zero. All 
other eigenvalues have negative real part. 

2. For an off-diagonal block {a — PojC-Po.fe} where there exists an intertwiner 
UPj — PiU with U e {PqHPq, PohaPo}' , there exists exactly one eigen- 
matrix to the eigenvalue zero. The eigenmatrix is Upj, where pj is the 
stationary eigenmatrix with support in PajH. All other eigenvalues have 
negative real part. 
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3. For an off-diagonal block {a = PojcxPo^k} where there exists an intertwiner 
U G {PqHPq — {EjPqj + EiPqi), PohaPo}' , there exists exactly one eigen- 
matrix to an eigenvalue on the axis of imaginary numbers. The eigenvalue 
is i{E£ — Ej) , the eigenmatrix is Upj, where pj is the stationary eigenmatrix 
with support in PojTi.. All other eigenvalues have negative real part. 

4. For an off-diagonal block {a = PojcrPo^k} where there is no intertwiner as 
in item\^ or\^ there exist only eigenvalues with negative real part. 

Proof. The existence and uniqueness of an eigenmatrix to the eigenvalue zero 
in cases ([T]), ([2]), and the nonexistence in case Q are stated and then proven 
in Theorem [13] and in Proposition 16; then Lemma 17 states that there are no 
generalized eigenspaces to this eigenvalue. 

It remains to examine the existence or non-existence of other eigenvalues on 
the imaginary axis. The method is the same as in the proof of Theorem 17 
"No circular paths" in |BNT08j . We switch between the Schrodinger and the 
Heisenberg picture. has the same spectrum as V. Assume the existence 
of an eigenvalue X = ir with eigenoperator F = PqiFPqj, where r G 



so 



exp [tT>'^)F = e^^^F, and use the Kadison inequality. No further details for the 



time evolution are needed to deduce equation (30). We refer to |BNT08j for 
description of how to conclude that V := F/ ^y\\F'^F\\ is a local isometry between 
PajH and PQ/Ti or a local unitary if j = i: 

V^-V = Poj, V-V^ = P 



0,i- 



(30) 



Now we use the evolution equation (12) for V, multiply from the left by V\ and 

get 



V^V\V) = t{V^HV-Po,jH) + J2{y^hiVha-lPo,jhiha-lV^hihaV) = irPo 



(31) 

Since, by definition of the collecting basins, {Poj^-Po/} C {Po-f^-fo, -Po^a-Po}', we 
may define 



Pojha — haP( 



0,j J 



V^KV, H = PojH, H = V^iiV, 



and write the trace of (31) as 



zTr[i^ - if] + J2mhiha] - lTT[hiha] - '^TilhiL]) = tr dim{Po,jH). (32) 



The Cauchy-Schwarz inequality, applied to inner product Tr[. ■ .] in the H.S. 
space of matrices, and the inequality between geometric and arithmetic mean 



imply that the real part of the l.h.s. of (32) is zero iff Va : ha = ha, i.e. 
V^haV = Pojha, haV = Vha- Usiug again (31 ) this implies H -V = V ■ {H + r). 



There are no such imaginary eigenvalues, no circular paths, in cases Tp|4 
They exist only in case|3| with U = V-\-V'<-\- @^^j ,^^i,Po^k- 



□ 
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3.6 Stationary states, collection of results 

A collection of results stated above in this Section gives now the proofs of the 
main theorems. 

Proof. For Theorem [2] on structuring of the Hilbert space: In Proposition 



14 on decomposition into basins the decomposition of TC into PoTi., defined as the 





lowest level of the cascade of decay, and its complement PqH is performed. This 



gives part (1) of the Theorem. 



The same Proposition 14 gives also the further splitting of Pq into basins. 
These are minimal collecting subspaces. If there exist stationary phase relations 
between two basins, and only then, can this splitting be varied, using another pair 
of basins which are related to the former two and among themselves by unitary 



transformations. This is stated and then proven in Proposition 16 on stationary 
phase relations. Their existence comes with a form of equivalence between basins, 
given by the unitary transformation, which has to commute with every Poh^Po 
and with PqHPq. 

There may exist a generalized form of equivalence between basins, where 
the unitary transformation commutes again with every P^h^Po but then with 
PqH Pq — EjPq j + EiPq i instead of PqHPq. Such a generalized form of equivalence 
implies again equivalence of the unique stationary states located in the basins; 
and it enables the occurrence of undamped oscillating phase relations. Collect- 
ing equivalent and generalized equivalent basins into subspaces Qo.fc^ gives the 
larger, unique, part of splitting Pq = J2k Qo,k- There are neither stationary nor 
undamped oscillating phase relations Qo^k'HQo/. This is stated and then proven 



in Theorem [18] on dephasing and non-dephasing. 

The recollection of basins can be seen as Qq^uH = C"(A;)(S)7ioo,fc- Uniqueness of 
the stationary states pk in the minimal collecting subspaces, everyone equivalent 



to ?ioo,fc! is stated in Theorem 13 So parts (2) (3) and (4) of the Theorem are 
proven. 



The splitting of -Pj^ i^^^ higher levels and basins is stated in Proposition 



16 



proving the last part of the Theorem. □ 

A subspace /C is the minimal support of a minimal stationary state iff it fulfills 
the conditions 



'Laziness" - the projector Pfc fulfills equation (22) 



'No creation of off-diagonal elements" - Pjc fulfills equation (24) 



• "Minimality" - /C contains no smaller subspace fulfilling the first two con- 
ditions. 

The characteristic equations have been found also by B. Kraus et al., |K08a[ 
IKOSb] . as determining "dark states", which are pure stationary states. 
A collection of results gives also the 
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Proof. For Theorem [3] on enclosures and blocks: The relations of projectors 
Qm G {'H, ha}' to the time evolution and to "enclosure" are analyzed in Lemma 
[7| and in the further discussions in Section 3.2 on "Conservation of projectors 
and enclosures". Because of the enclosure, all further decompositions can be 
performed for the restricted evolution acting on the set of density matrices p 
with support in Qm'H. □ 



4 Invariance (conservation) and symmetry 
4.1 Invariant observables 

The stationary states span a subspace of the HS-space of matrices, the space of 
eigenmatrices to the eigenvalue zero of the superoperator T>. Its adjoint, T>\ has 
the same set of eigenvalues, so there is a linear set of invariant operators. It 
has the same dimension as the set of stationary states, and contains only proper 



"eigenmatrices", in duality to Lemma 17 



Let us start a construction of an invariant observable with a nucleus, located 



in a collecting basin, say PQkti. The defining condition (23) "nothing goes out 



is equivalent to the dual condition "no observable comes in" 

Po,kX^\F)PQ^k = PQ,kT^\Po,kFPQ^k)Po,k- (33) 

Inside the basin is full dissipation with only one stationary state. There is there- 
fore only one invariant observable inside the basin (up to constant factors). It is 
the projector Pq ^. 

Po,kT'\Po,k)Po,k = Po,k. (34) 

The Heisenberg evolution goes backward. It lets, applied to the projector Po,^, the 
diagonal block Po,fc unchanged, it creates an extension into the decaying subspace 
P(j-, and also phase relations between this subspace and the basin where it started. 
It lets the whole block PqTC unchanged, since this collecting subspace is spanned 
by collecting basins, each one showing the "no observable comes in" condition 



(33). The evolved observable therefore stays HS-orthogonal to eventually existing 



undamped oscillating phase relations, and we can define 

Ao,fc := lim r*t(Po,fc) (35) 



t— >oo 



as an invariant observable. It is positive, since preserving positivity goes over 
from T* to T*^ by duality. 

For a maximal set of mutually orthogonal collecting basins one obtains 

5^Po,fc = Po ^ 5^Ao,fc = l. (36) 
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There may be a still larger set of linearly independent invariant observables. 
Let us represent subspaces with n equivalent collecting basins Pq/ allowing for 
stationary phase relations, but not undamped oscillating ones, as 



e 

For any n x n matrix M the observable M (g) 1 can be extended, analogously to 



the procedure (35), to an invariant observable, of these observables can be 



linearly independent, n of them are as constructed in (35). 

An algebraic analysis of the set of stationary states, of invariant observables 
and their relations to the generating operators relies on first cutting off the decay- 
ing subspace P^Ti,. The restricted time evolution, generated by V defined with 
Hamiltonian and transfer operators {PqHPq, Po^«-Po}; is identical to the full one 
for p with support on PqTY. For this restricted system, on PqH, the concepts 
of "basin" and "enclosure" are identical. Projectors onto basins are elements of 
the von Neumann algebra A/q = {PqHPq^ P^haPo}' ■, which is the set of invari- 
ant observables for the restricted time evolution. The set of extended invariant 
observables in the large system is thus related to the commuting von Neumann 
algebra of the restricted system; restricted to the non-decaying level PqH. 

3P: V\F) = F = Q ^ PoFPoeAfo (37) 

So, in systems without decay, the invariant observables do form an algebra. But, 



in systems with decay, the extension (35) does in general destroy this property. 
See examples in Section |5] 

4.2 Symmetries 

A dynamical symmetry is defined by existence of unitary or anti-unitary operators 
V, such that 

Vt, Vp T\VpV-^) = Vp{t)V-\ (38) 

This is equivalent to 

Vp V{VpV-^) = VV{p)V~\ (39) 

Sometimes an appearance of a dynamical symmetry is in connection with exis- 
tence of an algebra of conserved observables, a connection well known in Hamil- 
tonian dynamics. HUE {Hjha}' then T^{UpU'^) = Up{t)W. Another way how 
a symmetry may be guaranteed is, that VHV~^ = H and the set {VhaV~^} 
equals the set {ha}. But it may also be hidden, not immediately to be observed 
in the ha- See examples in Section [5j For continuous groups, and for T) with a fi- 
nite number of transfer operators ha-, the symmetry is necessarily not completely 
represented by invariance of the set of generators. 
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The maximal symmetry is invariance under all unitary and anti-unitary trans- 
formations. There is only one ray of generators compatible with this sym- 
metry, in the center of the cone of generators, see |BNT08j . It consists of 
{V\V{p) = X ■ {u — p)}, where u = 1/ dim(H) is the completely mixed state. 
There are several ways of choosing the set {ha} to form such a special T>. Exam- 
ples are again in Section |5| 

While, in Hamiltonian dynamics, the appearances of 

• dynamical symmetry 

• invariant (conserved) observables 

• algebra formed by the invariant observables 

must come together, these relations are not strict in irreversible dynamics. Here 
we observe cases of 

• dynamical symmetry without invariant observables, 

• Invariance of observables without a symmetry, 

• Invariant observables which do not form an algebra. 

For each of these cases we present examples. 

There remains a relation between a dynamical symmetry and the set of sta- 



tionary states. If is a symmetry operator as in (38), then the set of stationary 
states obeys the symmetry for stationarity 

{p\ stationary} = {V"pV"^^| p stationary}. (40) 



But this is, in general, a one-way relation. The symmetry for stationarity (40) 



may be valid, without (38) being true. This comes, again, because of the restric- 



tion of the set of invariant states to the subspace PoTi.. 



5 Examples 

The density matrices are representations of states in some basis of the usual type, 
employing a complete orthogonal set of basis vectors. Mostly we use V given as 
a sum of two simple generators, with transition operators {/i+,/i_}. Zeroes as 
matrix elements are represented with dots. Lower indices on C indicate the role 
of a subspace as a level or as a basin. Matrix elements of the density operator are 
denoted as Vij. Representations using a tensor product are in accordance to its 
use in part (3) of Theorem |2} "No invariant observable" means, precisely: Only 
the constants are invariant. 
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5.0.1 Dissipation 

Hilbert space H^C^ 

There is a unique stationary state. It is a; = 1/2, a fact which can not occur 
with simple generators. To represent the dynamics we use Pauh matrices, so 
h± = a± = ((Tx- + o'y)/2j and V{a±) = — o"±, T>{az) = —lo^- The dynamics 
is symmetric under rotation around the z-axis, under reflection cx^ —a^ and 
under complex conjugation C. The symmetry for stationarity is maximal, but 
there is no invariant observable. 



5.0.2 Decay, two collecting basins, no stationciry phase relations 

7i = = Co,i e Co,2 e Ci 

\ . . . \ /i_ = I . . 1 



Extremal stationary states are 1 © © and © 1 © 0. The symmetry for station- 
arity includes exchange of these extremal states and complex conjugation. The 
whole system is only one enclosure. Invariant observables are diagonal matrices 
with components (1,0,2/3), or (0,1,1/3), and their linear combinations. They 
do not form an algebra. No dynamical symmetry (but C). 



5.0.3 Decay, two collecting basins with stationary phase relations 

n = <c^ = Co,i © Co,2 © Ci 




Stationary states: Any 2x2 density matrix with support on Co,i®Co,2- Symmetry 
for stationarity: f/ © 1 with any unitary ?7, and C . No dynamical symmetry but 
C. Invariant observables: Linear combinations oi . . . A/^^ 




They do not form an algebra. 
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5.0.4 Decay, basins with dissipation inside 



/ ■ ±1 1 
1 ■ ■ 

V ■ ■ 1 



±1 



The restricted evolution on Cq is almost the same as in 5.0.1, only with twice 



the speed. One stationary state, a; © 0, no invariant observables. Symmetry for 
stationarity: all the unitary and anti unitary transformations acting on Cq. No 
invariant observables. 

5.0.5 Dephasing of two enclosures 

n = c^ 



(p2 



(p2 
^0,2 



/ ■ 1 ■ ■ \ / 



V- ■ - 

The restricted evolutions on each Cqj is the same as in 
off diagonal blocks, these evolutions are "out of phase", destroying every phase 
relation. Such an off diagonal block evolves according to 



5.0.1 



Acting on the 



d 
dt 



n,3 ^1,4 
^2,3 



^2,4 - ^1,3 
-^"2,3 



-ri,3 - r2,A 



Two extremal stationary states a; © and ©tu, two invariant observables, 1 © 
and © 1, very rich symmetry. 



5.0.6 Undamped oscillating phase relation 



H 



( 1 



V 



1 



/■ 1 


■ \ 




( ■ 








= 


1 






1 




V- ■ 


■ / 




V ■ 


1 ■/ 



The evolutions of states of the enclosures Cq • are again as in 



5.0.1 



Here they are 



in phase when acting on the off diagonal blocks, leaving a special phase relation 
undamped. But the Hamiltonian creates an oscillation. 



_^ f n,3 ^"1,4 

dt I r2,3 r2,4 



n,3 n,4 

r2,-i ^2,4 



+ 



^"2,4 - ri,3 -ri,4 

-^"2,3 ^"1,3 - ^"2,4 
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So 



1 1 n,3 ^^1,4 



~ e 



1 
1 



The equation for full dynamics can be written as 
d 



dt 



(M ® p) = -i[H, M] ® p + M ® Voflip), 



with Po,o as in 5.0.1 Stationary states are the same as in 5.0.5, They can be 



represented as tensor products of diagonal 2x2 matrices M with u. Symmetry of 
stationary states includes reflection, exchanging pi and p2- Invariant observables 
are given by all linear combinations of the two projectors onto Cq^. They form 
an algebra. 



5.0.7 Stationary phase relation 



The transition operators h± above in |5.0.6[ but V is given without the 

Hamiltonian. Very rich symmetry, much more symmetry operations than in 



5.0.5, There are more stationary states than in 5.0.5, they can be represented 
as tensor products M ® with any 2x2 matrix M > 0, Tr[M] = 1. Invariant 
observables are A ® 1; they form an algebra. 



5.0.8 Cascade of decay 

7f = = Co © Ci,i © Ci,2 © Ca 



/ ■ 1 ±1 • \ 



v 



■ / 



There are two lines of flow: |1,1)(1,1| ^ |0)(0| and |2)(2| ^|1,2)(1,2| ^ |0)(0|. 
Both lines have the same end. The differential equations are ri i = r2^2 + 2r3_3 
f2,2 = —''^2,25 ''^3,3 = — 2r3^3 + r44, 7^44 = — r44 for the diagonal matrix elements. 
The off-diagonals are just exponentially decaying, with one accompanying part 
of the flow, ri,3 = — 2ri^3 + r2,4 and the same for the adjoint. No invariant 
observables. 



5.0.9 Maximal symmetric evolution 

= C"; The evolution is 

p = dim(7Y) ■ {u 



P)- 



T> can be defined with H = 0, {hij = for some basis {|^)}. There is 

a unique stationary state, it is = 1/ dim(7Y). Symmetry under every unitary 
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Figure 1. Two cascades of decay. Basins are represented by rectangles, decaying as yellow, 
collecting as blue. Each flow between two basins is indicated by an arrow. 



and anti-unitary transformation holds. Another way to represent this evolution 
is to choose the Weyl operators as {hij} jBA08] . 

This is a special case of detailed balance at infinite temperature (see the remark 
around equation (2.15) in |G78j and references therein), which appears in all those 
cases, where {hl^} = {ha}. (Each ha either has a dual companion hp = h]^, or 
it is self-adjoint h]^ = ha-) One consequence of these symmetries of detailed 
balance is the invariance of the completely mixed state u = 1/ dim(?-^). Another, 
related, consequence is the self-adjointness of the superoperator V, ii H = 0. The 
examples 5.0.1 5.0.6 and 5.0.7 are cases of detailed balance. 



6 Perturbations 

6.1 Abstract considerations 

We consider generators V depending on a parameter A. The dependence on A of 
the transition operators is assumed as linear, which makes a quadratic dependence 
on A for Vx- We set Hx = Hq + XV + X'^W, and ha,x = ha + Xka- This gives 
Vx = Vq + X8 + with 

Po(p) = -?[i^,p] + $^{/^.p/^l-|(/^J./iaP + p/iL^a)}, (41) 

OL 

E{p) = -i[V,p] + Y,{hapki + kaphi 

a 

-\{hikaP + phika + kihaP + P^/l^)}, (42) 

Hp) = -imp\ + Y,{ko.pki-\{kikaP+ pkika)}. m 
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This formalism includes cases where the set of transition operators is enlarged. 
Formally, it is done by perturbing some = with 7^ 0. 



A general fact is that the defining equations for structuring, i.e. (22) for 



"laziness" and (24) for "nothing goes out", and for symmetry, i.e. (39), can 
abruptly turn to inequalities through an infinitesimal change A — >^ A + dX, but 
not the other way round: 

19 PROPOSITION. Consider the projection operators -P(A) onto subspaces, 
which are either enclosures, or collecting subspaces, or basins in the cascade of de- 
cay. Consider also the unitary and antiunitary operators V{X) which are symme- 



try operations as defined in Section 4-2. The functions A -P(A) and A 1-^ ^(A) 



are continuous functions, defined on closed sets 0/ A G M. 

Proof. To study "laziness" under the influence of perturbations ha ^a(A) 
define /(P, A) := J2a\\^a{X)P — Pha{X)P\\ for orthogonal projectors P. The 
projectors P can be expressed by a finite set of parameters, e.g. its matrix 
elements. The function /(P, A) is jointly continuous in P and A, so the set 
{(A, P)| /(P, A) = 0, P = P^ = P^}, which defines lazy subspaces via the implicit 
functions P(A), is closed. It is compact, when A is restricted to a compact 
interval. The implicit functions P(A) can be multi-valued; they give minimal 
lazy subspaces and direct sums of them. The domain of A G M for one P(A) is 
closed, the projector may disappear under infinitesimal changes of A. 

To study "collecting" subspaces, consider the functions /c(P, A) := \\P{iH{X) — 
^'l2a^aW^ce{X))P^\\, aud procccd in the same way as for laziness. To study 
symmetry, use fs{V,X) := \\VxiVpV-') - VD,{p)V-'\\. 

□ 

Infinitesimal changes of the transition operators and of the Hamiltonian can 
lead to mergers, may disturb an existing structure, or move it in the Hilbert 
space, but they can not create a new one. Moving basins in Hilbert space can 
not occur through mere addition of new transition operators, since the condition 



(22) for laziness involves each single ha- 
A list of things that can happen: 

• The number of zero eigenvalues may decrease (but it can not increase). 



Example: Disturb 5.0.2 with k± = ±|0, 1)(0, 2|. 



Stationary phase relations may turn to undamped oscillating ones, or they 
may become unstable. The inverse changes are not possible. Example: Dis- 



turb 5XL3 with the Hamiltonian = ±|0, 1)(0,1| or with k± = ±|0,1)(0,1|. 



Oscillating phase relations may disappear. Example: Disturb 5.0.6 with 



the h± of 5.0.5 



Enclosures may merge. Example : Disturb 5.0.5 with fc+ = |0)(1| (8> 1. 
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Collecting basins may merge. Example: The same as for the first item. 



Collecting basins may merge with decaying basins. Example: Disturb 5.0.8 
with fc± = ±|0)(1,2|. 

Rotations of basins in Hilbert space. Example: Consider the simple gener- 
ator with transition operator h = ax, disturb it with k = ay. 

Symmetries may disappear (but new symmetries can not emerge). Exam- 



ple: Below, in 6.2.1 the worked out perturbation of 5.0.1 



The general perturbation theory of linear operators, |Katoj . on finite di- 
mensional spaces states the analyticity of eigenvalues and eigenprojectors onto 
eigenspaces, with only algebraic singularities at some exceptional points. Now 
we are interested in the eigenvalue zero. Its multiplicity may decrease under per- 
turbation, but the eigenvalue zero has to remain, with multiplicity one at least. 
We want to follow this eigenvalue and the remaining eigenprojectors onto the 
stationary states. The constant function A i— *• is obviously analytic. So we have 
analyticity of the eigenprojector, from which we pick out the projectors onto 
those eigenmatrices which are states. Relying on the analyticity at A = for 

A ^ {Hy^,K,x} ^ Dx ^ {p{\)} 

we make an ansatz, expanding stationary states: 

oo 

P(A)=P + 5^AV.. (44) 

n=l 

We demand Tr[p(A)] = 1, which gives 

Tr[p] = 1, Tr[a„] = 0. (45) 

Expansion of the eigenvalue equation Vx{p{\)) = gives the starting condition 

Voip) = 0, (46) 

and a series of equations to determine the following contributions 

Vo{an) = -Sian-i) - ^(a„_2), (47) 

using cr_i = and ao = p. Trying to solve one of these equations one encounters 
two problems: Solving it requires that the r.h.s. is a matrix in rangefT'o]- Note 



that the range of all the operators involved in (47) consists of matrices with trace 



zero. This is part of solution to this first problem. The second problem: an 



inverse of Vq is not uniquely given, if (46) allows for more than one solution, i.e. 



if more than one stationary state exists. Choosing the right solution appears in 
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the iterated equations as solving there the first problem, deciding whether the 



r.h.s. of (47) is in the range of Vq. 

In order to proceed with calculations we define an inverse of Vq independently 
of the other super-operators: 

range[X'Q ^] := range[Po] = {'^1 which are invariant: Tr[Ajr] = 0}, (48) 

with using a maximal set of linearly independent invariant operators Ai. This 
makes the operator Vq^ unique: If there exists a, such that Vo^a) = r, then also 
T>Q{a + ajPj) = T, with HS-orthogonal eigenmatrices - stationary states and 
stationary phase relations - pj. Choose the Ai in such a way that Tr[y4jPj] = Sij, 
then 

TT[A{a + J2 "jPi)] = Tr[A*a] + = (49) 
j 

determines the uniquely. 

With this generally defined inverse of Vq one can invert mfh to 



cr„ = -Pq ^(^^(o"n-i) + ^(o-n-2)) + ^ftj-Pi, with ^ajTr[pj] = 0, (50) 

3 j 

and the task is now to find the right coefficients aj, so that the insertion of (T„ into 
the next iteration formulas gives matrices in the domain of definition of Vq^ - 
which is the range of T>o. Here we refrain from establishing a complete formalism, 
we represent several case studies instead, including special formulas. 

6.2 Case studies and examples 

6.2.1 Only one stationary state for the unperturbed system 

If there is only one p as eigenmatrix to the zero-eigenvalue of Vq, there are no 
more tasks to fulfill than to perform the calculations. It is not necessary to make 
a difference between systems with decay and those without. One may define 
superoperators Qn as 



and use them to calculate 



P(A)=P + VA"g„(p). (51) 



As an example to demonstrate the validity of this procedure we disturb |5.0.1 
with fc+ = (1 — crz)/2, leaving undisturbed. This gives a2n = and a2n+i = 
{—2)~^ax/2. The series converges for |A| < y/2 and can be summed up, giving 

P(A) = ^ + r:^^^-/2- (52) 
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6.2.2 Remaining enclosures, each one with only one stationary state 

Consider the case of projectors Qj commuting with Ha{X) and with every /iq(A). 



Inside each enclosure QjTi. the situation is the same as above, in 6.2.1 If there 
exist enclosures which allow pairwise stationary phase relations, collect them as 
C" (S> 7i. The evolution takes place only in Ti. and, if also the perturbation shows 



this symmetry, it can again be treated as in |6.2.1[ If it disturbs this symmetry, 
the phase relations vanish by dephasing. 

We remark that besides the families of stationary states Ylj PjPjW found in 
this way there could, formally, also be families with the Pj depending on A. 

6.2.3 Merging of enclosures through direct dissipation 

Consider H = Hi Q) H2, where each Hj = QjH is an enclosure supporting just 
one stationary state pj, and where no stationary phase relations exist. Perturb 



a (virtual) ha = with k, so the S in (42) is zero. The perturbing transition 



operator connects the enclosures via QikQ2 7^ 0, the other off-diagonal block of 
the matrix k may be zero or not. We demand moreover 

Qikp2 ^ 0. (53) 

The starting condition p = api + (1 — a)p2-, and the first of the conditions 



(47), which is J^{p) G range(r'o)) here to be expressed as Tr[QiJ^(p)] = 0, give. 



using Qi = 1 — Q2 and Ti[F{pi)] = 0, 

aTr[Qi.F(pi)] + (1 - a)Tr[gi^(p2)] = 
i:i[Qikp2k^Qi]-a{Tx[Q2kprk^Q2]+Tx[Qikp2k^Qi]) = 0. (54) 



That determines a G (0, 1] uniquely, since (g implies Ti[Qikp2k:^Qi] > 



To solve the following conditions of (47), 

^{(72n+2) = -J^(I'o ^(•^(o"2n))) + « " -^(pi " P2) G range(Po), 
we define the functional 

^ ^ TT[Q2kp,k^Q2] + TT[Qikp2k^Qi] ' 
The sequence an is now fixed as a2n+i = and 

(T2n+2 = -DQ^{T{a2n)) + (Pl - P2) " a[(X2n]- 

6.2.4 Dephasing perturbed by a Hamiltonian 

Consider again Ti = Tii ® 0.2, where each Hj = QjH is an enclosure support- 
ing just one stationary state pj, and where no stationary phase relations exist. 
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Perturb Vq with a Hamiltonian V, so the JF in (43) is zero. The perturbing V 



connects the enclosures via Q1VQ2 and Q2VQ1; we demand 

QiVp2 ^ 0. (55) 
For simphcity, assume QiVQi = and also Q2VQ2 = 0. 



The condition (55) implies, as can be checked by considering matrix elements 



involving eigenvectors of the pj , 

[V,pi-p2]j^0. 

Since Vq does not mix the matrix blocks QjaQk and does not annihilate 
phase relations, the range of Vq contains every QiaQ2, every and the 

inverse does exist for every off-diagonal block. This applies to S{p), where 
p = aoPi + (l — ao)p2- So the first order of perturbation theory gives no restriction 
on oq. We proceed with 



(Ti = -Vq \£ip2)) - oigDq ^{£{pi - P2)) + ai(pi 



P2J 



The condition £{ai) E range(r'o) does not involve ai but it determines the right 
value for ao, by demanding 

aoTr[gi^(I?o^(^(pi -P2)))] = -Tr[gi£(Po ^(^(pa)))], 
in case Tr[Qi ■ £(T)q^(£(pi — P2)))] 7^ 0. The abstract considerations of Sec- 



tion 6.1 imply that this factor has to be non-negative. In fact, the part due to 
— P2 may, with some tricky methods, be written as Tr[y4Qy4jj,], with = 

haT>Q^ {V P2) P2^^'^ Pp — Vq^ {V P2) P2^ Pphapl^'^ , whcrc Pp is the projector onto the 
range of p, and the ha are the transfer operators appearing in Vq. The part of 
the factor involving pi is zero if piV^ = 0, but in general it is also non-negative. 
We assume now, that at least one Aa 7^ 0, or, for short, just that the factor of 
ao, which appears also as factor for all the following a„, is not zero. 
The expansion of p(A) proceeds with 

cr„+i = 'DQ^{£{VQ^{£{an-i)))) - an'D^^{£{pi - P2)) + a„,+i(pi - P2)- 

The condition £{an+i) G range('Do) can be fulfilled by choosing = tt[o"n-i], 
with the functional 

-Ti[Q, ■ £{V,\£{V,\£{a)))))] 
i:v[Q,-£{V,\£{p,-p2)))] ■ 

6.2.5 Basins of a cascade which merge by dissipation 

Consider Ti = Tii® I-L2, where each Tij = QjH is spanned by the + 1 lowest 
levels ipjfi . . . ipj^N of a harmonic oscillator. The decay is performed through the 
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PerUorbation case 6.2.5 
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Figure 2. A perturbed cascade. Two independently decaying systems are perturbed in a way 
as is indicated by the broken pointed lines in red. 



annihilation operators aj acting as transition operators. This system is now 
perturbed by adding made of three new V. Two of them involve the creation 
operators Oj, the third one creates a dissipation between the two levels ipj^N- 
Finding the right constants a2n to do the expansion 

requires now an A^-fold iteration of applying the super-operator Q := —"Dq^ o JF. 
It takes steps of applying Q until the unperturbed ground states pj are lifted 
to the N'^^ level, and one more action of J-' to dissipate over to the other Tij. The 
procedure to get a2n+2 = a[c"2n] now involves the functional 

^ Tr[gi-(^og^+^)(a)] 
^ ^ Tr[gi-(^o^^)(pi-p2)]- 
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7 Summary 

Analysis of Lindblad generators is to be done on two levels. The lower level is 
the finite dimensional Hilbert space H spanned by the pure-state- vectors. The 
generators and their processes act at the upper level which contains the mixed 
states represented by density matrices. In this paper we have established an 
interplay of these two levels, where every process corresponds to a structure of 
Tt, a decomposition into mutual orthogonal subspaces. These subspaces, we call 
them basins, generalize the notion of energy levels which appears in Hamiltonian 
dynamics. Their dimensions may, in general, be any number between 1 and the 
dimension of the whole Hilbert space. They are mutually orthogonal, in spite of 
the possible non-hermiticity of the super-operators which represent the Lindblad 
generators and of the transfer operators acting on H.. 

The interplay between process and Hilbert space structure elucidates the char- 
acterization of the process, and brings about a way to structure the process itself, 
according to its action on blocks of density matrices. Decay corresponds to cer- 
tain subspaces which do not carry any stationary state. Minimal stationary states 
are supported by minimal collecting basins, inside of which Dissipation occurs. 
The phase relations between two basins, off-diagonal blocks of density matrices, 
either show Dephasing, or, in course of the process, go over to a special phase 
relation, which is either stationary or oscillating. 

Stationary and oscillating phase relations appear together with special dy- 
namical symmetries. The appearances of symmetries and of invariant obscrvables 
may show peculiar effects. In cases of Decay the invariant observables need not 
form an algebra. Dynamical symmetries are not the same as the symmetries 
appearing for the set of stationary states, in general. Perturbation of a process 
may lead to a merging of basins, which complicates the perturbative calculations 
of stationary states. 

Establishing the structure connected with a process and its Lindblad generator 
is, we think, a helpful tool for deeper investigations, probably indispensable. 



8 Appendix 

A block matrix form of operators to characterize subspaces 

The subspace which is to be characterized, is represented with the projector 
P = ^ Q Q ^ . It carries the density matrices P ~ q 
The subspace PTi. has the property of being 
Lazy ^h.^(^ 

(H G \ 
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an Enclosure -v^ H = ^ ^ ) ' ~ D 

See also |BNT08] . where formulas for calculating Vhip) in block matrix form 
are presented. 
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